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Abstract 

Brain is one of the most energy demanding organs in mammals, 
and its total metabolic rate scales with brain volume raised to a power 
of around 5/6. This value is significantly higher than the more com- 
mon exponent 3/4 relating whole body resting metabolism with body 
mass and several other physiological variables in animals and plants. 
This article investigates the reasons for brain allometric distinction 
on a level of its microvessels. Based on collected empirical data it 
is found that regional cerebral blood flow CBF across gray matter 
scales with cortical volume V as CBF ~ y^^/^, brain capillary diam- 
eter increases as V^^^"^, and density of capillary length decreases as 
y-i/6 predicted that velocity of capillary blood is almost in- 

variant (~ V"), capillary transit time scales as V^^^, capillary length 
increases as and capillary number as V'^^^~^, where e is typ- 

ically a small correction for medium and large brains, due to blood 
viscosity dependence on capillary radius. It is shown that the amount 
of capillary length and blood flow per cortical neuron are essentially 
conserved across mammals. These results indicate that geometry and 
dynamics of global neuro- vascular coupling have a proportionate char- 
acter. Moreover, cerebral metabolic, hemodynamic, and microvascular 
variables scale with allometric exponents that are simple multiples of 
1/6, rather than 1/4, which suggests that brain metabolism is more 
similar to the metabolism of aerobic than resting body. Relation of 
these findings to brain functional imaging studies involving the link 
between cerebral metabolism and blood flow is also discussed. 
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Introduction 

It is well established empirically that whole body metabolism of resting mammals 
scales with body volume (or mass) with an exponent close to 3/4, which is known as 
Kleiber's law [1, 2, 3, 4]. The same exponent or its simple derivatives govern the scalings 
of respiratory and cardiovascular systems in mammals and some other physiological 
parameters in animals and plants [2, 3, 5]. Because of its almost ubiquitous presence, 
the quarter power has often been described as a general law governing metabolism and 
blood circulation, and several formal models explaining its origin have been proposed 
that still cause controversy [6, 7, 8, 9]. However, as was found by the author [10], the 
brain metabolism at rest seems to follow another scaling rule. Total brain metabolic 
rate (both oxygen and glucose) scales with brain volume with an exponent ~ 0.85, or 
close to 5/6 [10]. Consequently, the volume-specific cerebral metabolism decreases with 
brain size with an exponent around —1/6, and this value is highly homogeneous across 
many structures of gray matter [10]. The origin of these cerebral exponents has never 
been explained, although it is interesting why brain metabolism scales different than 
metabolism of other systems. 

The brain, similar to other organs, uses capillaries for delivery of metabolic nutri- 
ents (oxygen, glucose, etc.) to its cells [11]. Moreover, numerical density of cerebral 
capillaries is strongly correlated with brain hemodynamics and metabolism [12, 13]. 
However, the cerebral microvascular network differs from other non-cerebral networks 
in two important ways. First, in the brain there exists a unique physical border, called 
the brain-blood barrier, which severely restricts infiux of undesired molecules and ions 
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to the brain tissue. Second, cerebral capillaries exhibit a large degree of physical plastic- 
ity, manifested in easy adaptation to abnormal physiological conditions. For instance, 
during ischemia (insufficient amount of oxygen in the brain) capillaries can substantially 
modify their diameter to increase blood flow and hence oxygen influx [14, 15, 16]. These 
two factors, i.e. structural differences and plasticity of microvessels, can in principle 
modify brain metabolism in such a way to yield different scaling rules in comparison 
to e.g. lungs or muscles. Another, related factor that may account for the uncommon 
brain metabolic scaling is the fact that brain is one of the most energy expensive organs 
in the body [10, 17]. This is usually attributed to the neurons with their extended axons 
and dendrites, which utilize relatively large amounts of glucose and ATP for synaptic 
communication [18, 19]. 

The main purpose of this paper is to determine scaling laws for blood flow and 
geometry of capillaries in the brain of mammals. Are they different from those found 
or predicted for cardiovascular and respiratory systems? If so, do these differences 
account for brain metabolic allometry? How the scalings of blood flow and capillary 
dimensions relate to the scalings of neural characteristics, such as neural density and 
axon (or dendrite) length? This study might have implications for expanding of our 
understanding of mammalian brain evolution, in particular the relationship between 
brain wiring, metabolism, and its underlying microvasculature [10, 20, 21]. The results 
can also be relevant for research involving the microvascular basis of brain functional 
imaging studies, which use relationships between blood flow and metabolism to decipher 
regional neural activities [22, 23]. 
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Results 

The data for brain circulatory system were collected from different sources (see 
Materials and Methods) . They cover several mammals spanning 3-4 orders of magnitude 
in brain volume, from mouse to human. 

1. Empirical scaling data. 

Cerebral blood flow CBF in different parts of mammalian gray matter decreases 
systematically with gray matter volume, both in the cortical and subcortical regions 
(Fig. 1). In the cerebral cortex, the scaling exponent for regional CBF varies from 
—0.13 for the visual cortex (Fig. lA), —0.15 for the parietal cortex (Fig. IB), —0.17 for 
the frontal cortex (Fig. IC), to —0.19 for the temporal cortex (Fig. ID). The average 
cortical exponent is —0.16 ±0.02. In the subcortical regions, the CBF scahng exponent 
is —0.14 for hippocampus (Fig. 2A), —0.17 for thalamus (Fig. 2B), and —0.18 for 
cerebellum (Fig. 2C). The average subcortical exponent is identical with the cortical 
one, i.e., — 0.16±0.02, and both of them are close to —1/6. It is interesting to note that 
almost all of the cortical areas (except temporal cortex) have scaling exponents whose 
95% confidence intervals do not include a quarter power exponent —1/4. 

The microvessel system delivering energy to the brain consists of capillaries. The 
capillary diameter increases very weakly but significantly with brain size, with an ex- 
ponent of 0.08 (Fig. 3A). On the contrary, the volume-density of capillary length 
decreases with brain size raised to a power of —0.16 (Fig. 3B). Thus, the cerebral 
capillary network becomes sparser as brain size increases. Despite this, the fraction 
of gray matter volume taken by capillaries is approximately independent of brain size 
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(Fig. 3C). Another vascular characteristic, the arterial partial oxygen pressure, is also 
roughly invariant with respect to brain volume (Fig. 3D). 

A degree of neurovascular coupling can be characterized by geometric relationships 
between densities of capillaries and neurons. Scaling of the density of neuron number in 
the cortical gray matter is not uniform across mammals [24, 25, 26]. In fact, the scaling 
exponent depends to some extent on mammalian order and the animal sample used [26] . 
For the sample of mammals used in this study, it is found that cortical neuron density 
decreases with cortical gray matter volume with an exponent of —0.13 (Fig. 4A). This 
exponent is close to the exponent for the scaling of capillary length density, which is 
—0.16 (Fig. 3B). Consistent with that, the ratio of cortical capillary length density to 
neuron density across mammals is approximately constant and independent of brain 
size (Fig. 4B). Typically, there is about 10 iim of capillaries per cortical neuron. The 
scaling dependence between the two densities yields an exponent close to unity (Fig. 
4C), which shows a proportionality relation between them. 

Cerebral blood flow CBF scales with brain volume the same way as does capillary 
length density (Figs. 1,2,3B), and thus, CBF should also be related to neural density. 
Indeed, in the cerebral cortex the ratio of the average CBF to cortical neural density 
is independent of brain scale (Fig. 5). This means that the average amount of cortical 
blood flow per neuron is invariant among mammals, and about (1.45±0.4)10~^ mL/min. 
Taken together, the findings in Figs. 4 and 5 suggest a tight global correlation between 
neurons and their energy supporting microvascular network. 
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2. Theoretical scaling rules for cerebral capillaries. 

Below I derive theoretical predictions for the allometry of brain capillary character- 
istics, such as: capillary length and radius, capillary number, blood velocity, and time 
taken by blood to travel through a capillary. I also find relationships connecting cere- 
bral metabolic rate and blood flow with neuron density. The following assumptions are 
made in the analysis: (i) Oxygen consumption rate in gray matter CMR02 scales with 
cortical gray matter volume V as V~^^^, in accordance with Ref. [10]; (ii) Capillary 
volume fraction, fc = ttNcLcRI/V, is invariant with respect to V, which follows from 
the empirical results in Fig. 3C. The symbol Nc denotes total capillary number in the 
gray matter, is the length of a single capillary segment, and Rc is its radius; (iii) 
Driving blood pressure Ap^ through capillaries is independent of brain size, which is 
consistent with a known fact that arterial blood pressure (both systolic and diastolic) of 
resting mammals is independent of body size [27, 28, 29] ; (iv) Partial oxygen pressure 
P02 in capillaries is also invariant, which is consistent with the empirical data in Fig. 
3D on the invariance of arterial oxygen pressure; (v) Cerebral blood flow CBF is pro- 
portional to oxygen consumption rate CMR02, due to adaptation of capillary diameters 
to oxygen demand. 

The cerebral metabolic rate of oxygen consumption CMR02, according to the modi- 
fled Krogh model [11, 14] , is proportional to the product of oxygen flux through capillary 
wall and the tissue-capillary gradient of oxygen pressure Apo2, i-e. 




(1) 
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where D is the oxygen diffusion constant in the brain. The dependence of CMR02 
on capillary radius in this model has mainly a logarithmic character, and hence it is 
neglected as weak. Since oxygen pressure in the brain tissue is very low [30] , the pressure 
gradient Apo2 is essentially equal to the capillary oxygen pressure po2- Consequently, 
the formula for CMR02 simphfies to CMR02 ~ PcP02, where pc is the density of capillary 
length pe = KLJV. 

From the assumptions (i) and (iv) wc obtain that capillary length density Pc ~ V~^^^. 
Additionally, from (ii) we have ~ fc/Pc ~ V^/V~^/^ ~ V^^^, implying that capillary 
radius (or diameter) Rc scales as V^^^'^. Consequently capillary diameter does not 
increase much with brain magnitude. As an example, a predicted capillary diameter for 
elephant with its cortical gray matter volume 1379 cm^ [31] is 7.2 /im, which does not 
differ much from those of rat (4.1 //m [15, 32]) or human (6.4 //m [33, 34]), who have 
corresponding volumes 3450 and 2.4 times smaller. 

The blood flow Qc through a capillary is governed by a modified Poiseuille's law in 
which blood viscosity depends on capillary radius [35] : 



where Ap^ is the axial driving blood pressure along a capillary of length Lc, and r]ef{Rc) 
is the capillary radius dependent effective blood viscosity. The latter dependence has a 
nonmonotonic character, i.e. for small diameters the viscosity f]ef{Rc) initially decreases 
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with increasing Rc, reaching a minimum at diameters about 5 — 7 fim. For 2Rc > 10 
/im the blood viscosity rjef slowly increases with Rc approaching its bulk value for 
diameters ~ 500 /im. This phenomenon is known as the Fahraeus-Lindqvist effect [36] . 
In general, blood viscosity in narrow microvessels depends on microvessel thickness 
because red blood cells tend to deform and place near the center of capillary leaving 
a cell-free layer near the wall [35, 37]. These two regions have significantly different 
viscosities, with the cell-free layer having essentially plasma viscosity rjp, which is much 
smaller than the bulk (or center region) viscosity r]c- The formula relating the effective 
blood viscosity rjef with capillary radius and both viscosities rjp and r)c is given by [35]: 



where = r]p/r)c <^ 1, and w is the thickness of cell- free layer. 

For capillary radiuses relevant for the brain, i.e. 1.5 /im < Rc < 3.5 /im (see Suppl. 
Table S2), the ratio w/Rc increases with increasing R^ which causes a decline in the 
effective blood viscosity down to its minimal value at i?c = 3 — 3.5 /zm (Table 1). Using 
the data in Table 1 taken from [35], we can approximate the denominator in Eq. (3) 
for this range of radiuses by a simple, explicit function of Rc- The best fit is achieved 
with a logarithmic function, i.e. 1 — (1 — — w/Rc)^ ~ 0.85[\n(Rc/ Rq)]'^^^ , where 
Ro — 1.2 /im (Table 1). As a result, the effective blood viscosity takes a simple form: 



Vef{Rc) 



l-{l-fj,){l-w/Rc) 



4' 



(3) 



ricfiRc) = l.lSr|p{ln{Rc/Ro)^ 



{^) 
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Cerebral blood flow CBF in the brain gray matter is deflned as CBF = Q/V^ where 
Q — NcQc is the total capillary blood flow through all Nc capillaries. Thus CBF is 
given by 



CBF 



VefiRc) LI ' 



(5) 



or 



CBF 



^PcPcR 



4 r 



ln(^ 
Ro 



2/3 



(6) 



We can rewrite the logarithm present in Eq. (6), in an equivalent form, as a power 
function [Rc/Ro)'^ with a variable exponent 7 given by (see Appendix SI in the Supp. 
Infor.): 



7=0 



2 1n(ln(i?,/ig„)) 
3 HRjRo) ' 



(7) 



so that CBF becomes 



LI 



(8) 



The exponent 7 in this equation can be viewed as a correction due to non-constant 
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blood viscosity (Fahraeus-Lindqvist effect [36]). Tfie dependence of 7 on the capillary 
diameter is shown in Table 2. Because in general 7 < 0, its presence in Eq. (8) 
reduces the power of Rc- However, this effect is weak for medium and large brains as 
I7I <S 1. Even for a small rat brain the relative influence of 7 is rather weak, since 
I7I/4 a; 0.19. In contrast, for very small brains, such as mouse, the effect caused by 7 
is strong (Table 2), which reflects a sharp increase in the effective blood viscosity for 
the smallest capillaries [35, 37]. 

Now we are in a position to derive scaling rules for the capillary length segment 
Lc, capillary blood velocity Uc, and the number of capillaries Nc- Prom Eq. (8), using 
the assumptions (i), (iii), and (v), we obtain ~ pcRf^'^ / CMRo2, which implies that 
Lc ~ (viscosity of blood plasma is presumably independent of brain scale [38]). 

Consequently, Lc ~ \/i/6+7/24^ j g capillary length should weakly increase with brain 
size. Although there are no reliable data on Lc, we can compare our prediction with the 
measured intercapillary distances, which generally should be positively correlated with 
Lc- Indeed, the mean intercapillary distance in gray matter increases with increasing 
brain volume, and is 17 — 24 /im in rat [39], 24 /im in cat [40], and 58 /im in human 
[33]. 

Average velocity Uc of blood flow in brain capillaries is given by Uc — Qc/{t^R^)- 
Using the expressions for Qc and r^e/, we get Uc ~ {Apf.Rl'^'^) / Lc- Since Lf. ~ R^'^'^^'^ 
above, and using the assumption (iii), we obtain Uc ~ RJ^'^ ~ ■\/7/24_ xhus, capillary 
blood velocity is almost independent of brain size for medium and large brains, as then 
7 — > (Table 2). For very small brains, instead, there might be a weak dependence. A 
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related quantity, the blood transit time Tc through a capillary, defined as Tc — Lc/uci 
scales as Tc ~ V^^^^ regardless of the brain magnitude. This indicates that Tc and 
CBF are inversely related across different species, Tc ~ CBF~^, because of their scahng 
properties. 

We can find the scaling relation for the total number of capillaries N^. from the 
volume-density of capillary length p^. We obtain = pcV/L^ ~ y-V6y/yi/6+7/24 ^ 
y2/3-7/24^ i.e. the exponent for Nc is close to 2/3 for not too small brains. As an 
example, the number of capillary segments in the human cortical gray matter should 
be 123 times greater than that in the rat (cortical volumes of both hemispheres in rat 
and human are 0.42 cm^ [24] and 572.0 cm^ [41], respectively). 

As was shown above, CMR02 must be proportional to the volume density of capillary 
length pc (Eq. 1). On the other hand, the empirical results in Fig. 4 indicate that pc 
is roughly proportional to neuron density pn- Thus, we have approximately CMR02 
~ Pn across different mammals. This implies that oxygen metabolic energy per neuron 
in the gray matter should be approximately independent of brain size. Exactly the 
same conclusion was reached before in a study by Herculano-Houzel [26], based on 
independent data analysis. Moreover, since cortical CMR02 and CBF scale the same 
way against brain size, we also have CBF ~ p„, which is confirmed by the results in 
Fig. 5. In other words, both cerebral metabolic rate and blood flow per neuron are 
scale invariant. 
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Discussion 



1. General discussion. 

The summary of the scahng resuhs is presented in Tabic 3. Some of these allomet- 
ric relations are directly derived from the experimental data (CBF, Rc, Pc, fc, Pc/ Pn, 
CBF/p„), and others are theoretically deduced {Nc, Lc, Uc, Tc). The interesting result 
is that cerebral blood flow CBF in gray matter scales with cortical gray matter volume 
raised to a power of —0.16. The similar exponent governs the allometry of cortical 
metabohc rate CMR [10], which indicates that brain metabolism and blood flow are 
roughly linearly proportional across different mammals. This conclusion is compatible 
with several published studies that have shown the proportionality of CMR and CBF 
on a level of a single animal (rat, human) across different brain regions [12, 42]. 

The coupling between CMR and CBF manifests itself also in their relation to the 
number of neurons. In this respect, the present study extends the recent result of 
Herculano-Houzel [26] about the constancy of metabolic energy per neuron in the brains 
of mammals, by showing that also cerebral blood flow and capillary length per neuron 
are essentially conserved across species. There are approximately 10 /im of capillar- 
ies and 1.45 • 10~^ mL/min of blood flow per cortical neuron (Figs. 4 and 5; Supp. 
Tables S2 and S3). This finding suggests that not only brain metabolism but also its 
hemodynamics and microvascularization are evolutionarily constrained by the number 
of neurons. This mutual coupling might be a result of optimization in the design of 
cerebral energy expenditure and blood circulation. 

It should be underlined that both CBF and CMR scale with brain volume with the 
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exponent about —1/6, which is significantly different from the exponent —1/4 relating 
whole body resting specific metabolism with body volume [1, 2, 3]. Instead, the cerebral 
exponent —1/6 is closer to an exponent —0.12 ± 0.02 characterizing maximal body 
specific metabolic rate and specific cardiac output in strenuous exercise [43, 44]. In 
this sense, the brain metabolism and its hemodynamics resemble more the metabolism 
and circulation of exercised muscles than other resting organs, which is in line with the 
empirical evidence that brain is an energy expensive organ [10, 17, 18]. This may also 
suggest that there exists a common plan for the design of micro circulatory system in 
different parts of the mammalian body that uses the same optimization principles [45] . 

The results of this study show that as brain increases in size its capillary network 
becomes less dense, i.e. the densities of both capillary number and length decrease, 
respectively as Nc/V ~ \/-V3-7/24 ^ y-i/a (Table 3). Contrary to that, the 

capillary dimensions increase weakly with brain volume, their radius as Rc ~ T/^/^^ and 
their length segment as ~ which are sufficient to make the fraction 

of capillary volume in the gray matter to be scale invariant (Table 3). The correction 
7/24 appearing in the scaling exponents for Nc/V and Lc reflects the fact that blood 
viscosity depends on capillary radius (Fahraeus-Lindqvist effect [36]). This correction 
is however small for sufficiently large brains, generally for brains larger or equal to that 
of rat, for which typical values of 7/24 are in the range from —0.032 to —0.0004 (Table 
2). On the contrary, for brains of mouse size or smaller, this correction is substantial, 
about —0.15, which implies that for very small brains Lc is essentially constant. 

Despite the changes in the geometry of microvessels, the velocity of capillary blood 
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Uc is almost scale invariant for not too small brains (exponent 7/24 f» 0; Table 3). This 
prediction agrees with direct measurements of velocity in the brains of mouse, rat, and 
cat, which does not seem to change much, i.e. it is in the range 1.5 — 2.2 mm/sec 
[40, 46]. Consequently the transit time Tc through a capillary increases with brain size 
as Tc ~ V^/^^ i.e. the scaling exponent is again 1/6. Another variable that seems to be 
independent of brain scale is partial oxygen pressure in cerebral capillaries (Table 3), 
which is consistent with the empirical findings in Fig. 3D on the invariance of oxygen 
pressure in arteries, as the two circulatory systems are mutually interconnected. 

2. Capillary scaling in cerebral and non-cerebral tissue. 

The above scaling results for the brain can be compared with available analogous 
scaling rules for pulmonary, cardiovascular, and muscle systems. For these systems, it 
was proposed (no direct measurements) that partial oxygen pressure in capillaries should 
decline weakly with whole body volume (or organ volume as lung and heart volumes, 
yiung-i Vheart, scalc isomctrically with body volume [2]) with an exponent around —1/12, 
to account for the whole body specific metabolic exponent —1/4 [47, 48]. In the resting 
pulmonary system, the capillary radius as well as the density of capillary length scale the 
same way as they do in the brain, i.e., with the exponents 1/12 and —1/6, respectively, 
against system's volume [49]. Also, the capillary blood velocity in cerebral and non- 
cerebral tissues scale similarly, at least for not too small volumes, i.e. both are scale 
invariant [2, 3] (Table 3). However, the number of capillaries and capillary length seem 
to scale slightly different in the resting lungs, i.e. Nc ~ Viuiig and Lc ~ Viung [47], 
although the difference can be very mild. For the resting heart, it was predicted (again. 
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no direct measurements) that ~ yhiarti blood transit time through a capillary 
r ~ yhiart i-^- the exponents are multiples of a quarter power and are slightly larger 
than those for the brain (Table 3). Interestingly, for muscles and lungs in mammals 
exercising at their aerobic maxima, the blood transit time scales against body mass 
with an exponent close to 1/6 [50], which is the same as in the brain (Table 3). This 
again suggests that brain metabolism is similar to the metabolism of other maximally 
exercised organs. Overall, the small differences in the capillary characteristics among 
cerebral and non-cerebral resting tissues might account for the observed differences in 
the allometries of brain metabolism and whole body resting metabolism. In particular, 
the prevailing exponent 1/6 found in this study for brain capillaries, instead of 1/4, 
seems to be a direct cause for the distinctive brain metabolic scaling. 

3. Brain microvascular network vs. neural network. 

The interesting question from an evolutionary perspective is how the allometric 
scalings for brain capillary dimensions relate to the allometry of neural characteristics. 
The neural density pn (number of cortical neurons per cortical gray matter volume 
V) scales with cortical volume with a similar exponent as does the density of capillary 
length pc (Fig. 4A). Thus, as a coarse-grained global description we have approximately 
Pn ~ Pc (Fig- 4B,C), or ~ N^Lc- The latter relation means that the total number 
of neurons is roughly proportional to the total length of capillaries, or equivalently, 
that capillary length per cortical neuron is conserved across different mammals. This 
cross-species conclusion is also in agreement with the experimental data for a single 
species. In particular, for mouse cerebral cortex it was found that densities of neural 
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number and microvessel length are correlated globally across cortical areas (but not 
locally within a single column) [51]. Moreover, since axons and dendrites occupy a 
constant fraction of cortical gray matter volume (roughly 1/3 each; [52, 53]), we have 
Nnld? ~ V , where I and d are respectively axon (or dendrite) length per neuron and 
diameter. Furthermore, because the average axon diameter d (unmyelinated) in the 
cortical gray matter is approximately invariant against the change of brain scale [52, 54], 
we obtain the following chain of proportionahties: I ~ ~ pj^ ~ V^/^ ~ L", where 
the exponent a = 1/(1 + 7/4). For medium and large brains, a ~ 1, implying a nearly 
proportional dependence of axonal and dendritic lengths on capillary segment length. 
For very small brains (roughly below the volume of rat brain), a can be substantially 
greater than 1, suggesting a non- linear dependence between capillary and neural sizes. 

Given that the main exchange of oxygen between blood and brain takes place in 
the capillaries, these results suggest that metabolic needs of larger brains with greater 
but numerically sparser neurons must be matched by appropriately longer yet sparser 
capillaries. This finding reflects a rough, global relationship, which might or might not 
be related to the fact that during development neural and microvessel wirings share 
mutual mechanisms [20, 55]. At the cortical microscale, however, things could be more 
complicated, and a neuro-vascular correlation might be weaker, as both systems are 
highly plastic even in the adult brain (e.g. [56]). Regardless of its nature and precise 
dependence, the neuro-vascular coupling might be important for optimization of neural 
wiring [53, 57, 58]. In fact, neural connectivity in the cerebral cortex is very low, and 
it decreases with brain size [58, 59], similar to the density of capillary length (Fig. 3B, 
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Table 3). To make the neural connectivity denser, it would require longer axons and 
consequently longer capillaries. That may in turn increase excessively brain volume 
and its energy consumption, i.e. the costs of brain maintenance. As a result, the 
metabolic cost of having more neural connections and synapses for storing memories 
might outweigh its functional benefit. 

The brain metabolism is obviously strictly related to neural activities. In general, 
higher neural firing rates imply more cerebral energy consumed [18, 19]. It was esti- 
mated, based on a theoretical formula relating CMR with firing rate, that the latter 
should decline with brain size with an exponent around —0.15 [19]. This implies that 
neurons in larger brain are on average less active than neurons in smaller brains. Such 
sparse neural representations may be advantageous in terms of saving the metabolic en- 
ergy [18, 60, 61]. At the same time, what may be related, neural activity is distributed 
in such a way that both the average energy per neuron and the average blood flow per 
neuron are approximately invariant with respect to brain size (Fig. 5; Table 3, [26]). 
Additionally, average firing rate should be inversely proportional to the average blood 
transit time Tc through a capillary, because both of them scale reversely with brain size 
(Table 3). Thus, it appears that global timing in neural activities should be correlated 
with the timing of cerebral blood flow. These general considerations suggest that apart 
from structural neuro- vascular coupling there is probably also a significant dynamic 
coupling. This conclusion is qualitatively compatible with experimental observations 
in which enhanced neural activity is invariably accompanied by increase in local blood 
flow [62]. 
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4. Relationship to brain functional imaging. 

The interdependencies between brain metabolism, blood flow, and capillary param- 
eters can have practical meaning. Currently existing techniques for non-invasive visu- 
ahzation of brain function, such as PET or fMRI, are associated with measurements of 
blood flow CBF and oxygen consumption CMRo2- It turns out that during stimulation 
of a specific brain region, CBF increases often, but not always, far more than CMR02 
[63]. However, both of them increase only by a small fraction in relation to the back- 
ground activity, even for massive stimulation [62, 63]. This phenomenon was initially 
interpreted as an uncoupling between blood perfusion and oxidative metabolism [64]. 
Later, it was shown that this asymmetry between CBF and CMR02 can be explained 
in terms of mechanistic limitations on oxygen delivery to brain tissue through blood 
flow [65]. We can provide a related, but simpler explanation of these observations that 
involves physical limitations on the relative changes in capillary oxygen pressure and 
radius. 

During brain stimulation, both CBF and CMR02 change by (5CBF and 5CMRo2, 
which are according to Eqs. (1) and (8) related to modifications in capillary radius 
(from Rc to Rc + SRc), and changes in partial oxygen pressure {po2 ^ P02 + Spo2)- The 
density of perfused capillary length pc remains constant for normal neurophysiological 
conditions. Accordingly, a small fraction of blood flow change is 



5CBF 
CBF 



(4 + 7) 



SRc 



Rc 



(9) 
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and similarly, a small fractional change in the oxygen metabolic rate is: 



CMR02 



PO2 



(10) 



In general, oxygen pressure increases with increasing capillary radius, in response to 
increase in blood flow CBF. This relationship can have a complicated character. We 
simply assume that po2 ~ -R", where the unknown exponent a (a > 0) contains all the 
non-linear effects, however complicated they are. Thus, a small fractional change in 
oxygen pressure can be written as 5po2/po2 ~ aSRc/Rc- As a result, we obtain 



If partial oxygen pressure po2 depends on capillary radius linearly or sublinearly, i.e., 
if a < 1, then the fractional increase in oxygen metabolism is signiflcantly smaller 
than a corresponding increase in cerebral blood flow. This case corresponds to the 
experimental reports showing that this ratio is <^ 1, for example, in the visual cortex 
(~ 0.1) [66] and in the sensory cortex (~ 0.2 — 0.4) [64, 67]. If, in turn, po2 depends 
on Rc superlinearly, i.e. if a > 1, then the coefficient a/(4 + 7) in Eq. (10) can be of 
the order of unity. Such cases have been also reported experimentally during cognitive 
activities [42] or anesthesia [68, 69]. 



CMR02 



a 6CBF 
(4 + 7) CBF ■ 



(11) 
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Materials and Methods 

The ethics statement does not apply to this study. CBF data were collected from 
different sources: for mouse [70], rat [71], rabbit [72], cynomolgus monkey [73], rhesus 
monkey [74] , pig [75] , and human [76] . Cerebral capillary characteristics were obtained 
from several sources: for mouse [14, 51], rat [15, 77, 32], cat [40, 78], dog [79], rhesus 
monkey [80], and human [33, 34]. Data for calculating neuron densities were taken 
from [24, 25, 41, 52, 81, 82]. Cortical volume data (for 2 hemispheres) are taken from 
[52, 41, 81]. Their values are: mouse 0.12 cm^, rat 0.42 cm^, rabbit 4.0 cm^, cat 14.0 
cm^, cynomolgus monkey 21.0 cm^, dog 35.0 cm^, rhesus monkey 42.9 cm^, pig 45.0 cm^, 
human 571.8 cm^. All the numerical data are provided in the Supporting Information 
(Tables SI, S2, and S3). 



21 



Supporting Information 

Appendix SI 
Table SI 

Regional cerebral blood flow CBF in mammals. 
Table S2 

Cerebral capillary and neural characteristics in mammals. 
Table S3 

Arterial partial oxygen pressure and average cortical CBF per neuron. 
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Figure Captions 
Fig. 1 

Scaling of cerebral blood flow CBF in the cortical gray matter. (A) Visual cortex: 
y = -0.127X + 1.98 {R^ = 0.93, p < 0.001). 95% confidence interval for the slope CI=(- 
0.168,-0.086). (B) Parietal cortex: y = -0.150a; + 2.00 {R^ = 0.89, p = 0.005), slope 
CI=(-0.222,-0.078). (C) Frontal cortex: y = -0.170x + 2.00 {R^ = 0.89, p = 0.002), 
slope CI=(-0.239,-0.100). (D) Temporal cortex: y = -0.191x + 2.09 (i?^ = 0.89, 
p = 0.005), slope CI=(-0.286,-0.096). 

Fig. 2 

Scaling of cerebral blood flow CBF in the subcortical gray matter. (A) Hippocampus: 
y = -0.135X + 1.89 {R^ = 0.90, p = 0.049), slope CI=(-0.271,0.000). (B) Thalamus: 
y = -0.167a; + 2.02 {R^ = 0.83, p = 0.011), slope Cl=(-0.272,-0.062). (C) Cerebellum: 
y = -0.177a; + 2.03 {R^ = 0.88, p = 0.002), slope CI=(-0.252,-0.102). 

Fig. 3 

Scaling of brain capillary characteristics against brain size. (A) Capillary diameter 
scales against cortical gray matter volume with the exponent 0.075 {y — 0.075a; + 0.60, 
i?^ = 0.87, p = 0.007), exponent CI=(0.034,0.117). (B) Volume density of capillary 
length scales with the exponent -0.16 {y = -0.162a; + 2.86, R^ = 0.79, p = 0.044), 
exponent CI=(-0. 316,-0. 008). (C) Fraction of capillary volume fc in gray matter is 
essentially independent of brain size {y = 0.029a; — 1.83, R^ = 0.07, p = 0.662), the same 
as (D) the arterial partial oxygen pressure {y — 0.012a; -|- 1.96, R"^ — 0.09, p — 0.472). 
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Fig. 4 

Neuron density versus capillary length density in the cerebral cortex. (A) Across our 
sample of mammals, the cortical neuron number density p„ scales against cortical vol- 
ume with the exponent -0.13 {y = -0.128x + 1.86, = 0.87, p = 0.022), exponent 
CI=(-0. 221,-0. 036). (B) The ratio of the density of capillary length to the density 
of neurons p„ in the cortex does not correlate with brain size (y = —O.OSAx + 0.99, 
i?2 = 0.09, p = 0.617), exponent CI=(-0.228,0.161). (C) The log-log dependence 
of the capillary length density pc on neuron density p„ gives the exponent of 1.05 
{y = imix + 0.88, = 0.63, p = 0.109). 

Fig. 5 

Invariance of cerebral blood flow per cortical neuron across mammals. The ratio of 
CBF to neuron density p„ in the cerebral cortex does not correlate significantly with 
brain size (log-log plot yields y = -0.033a; + 0.183, R"^ = 0.39, p = 0.261). The value 
of CBF for each species is the arithmetic mean of regional CBF across cerebral cortex. 
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Table 1: Parameters affecting the effective blood viscosity. 



Rc \pM 


w [//m] 


w/Rc 


l-{l-fi){l-w/Rc)^ 


0.85(ln(i?e/1.2))2/3 


1.5 


0.07 


0.05 


0.27 


0.31 


2.0 


0.30 


0.15 


0.54 


0.54 


2.5 


0.60 


0.24 


0.71 


0.69 


3.0 


0.90 


0.30 


0.79 


0.80 



Data for w and Rc were collected from [35]. The value of fj, was taken as 1/8. 
The last column represents values of the fitting function to the function in 
the fourth column. 



Table 2: Exponent 7 as a function of capillary diameter. 



Species 


mouse 


rat 


cat 


dog 


monkey 


human 


2Rc [//m] 


3.1 


4.1 


5.1 


4.5 


5.6 


6.4 


7 


-3.51 


-0.77 


-0.25 


-0.49 


-0.13 


-0.01 



Data for Rc were taken from Suppl. Inform. Table S2 (references therein). 



Table SI: Reg 


ional cere 


)ral 1)1()()(1 


flow CBF 


ill iiiaiiiiiialj- 






Species 


visual 


frontal 


temporal 


parietal 


hippocam. 


thalamus 


cerebell. 




cortex 


cortex 


cortex 


cortex 








Mouse ('^^ 


124.0 


165.0 


157.0 


166.0 


112.8 


195.5 


195.0 


Rat 


116.0 


124.0 


195.0 


113.0 


77.8 


116.0 


123.0 


Rabbit 


69.6 


67.2 




63.6 




60.6 


62.4 


Cynomolgus 
















monkey ^'^^ 


55.4 


48.7 


70.6 


56.8 




56.7 


57.2 


Rhesus monkey 


59.2 


45.0 


53.2 




53.9 




47.8 


Pig 


66.5 


63.3 


45.5 


62.8 


41.0 


48.5 


64.0 


Human 


43.0 


41.2 


44.6 


43.4 




47.5 


41.5 



CBF data are given in mL/(100g*min). 

References: ('^) [70]; [71]; ('^) [72]; ^ [73]; [74]; [75]; [76]. 
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Table 3: Summary of scalings for brain capillaries and hemodynamics against 
cortical gray matter volume V. 



Parameter 



Scaling rule 



Capillary radius, Rc 
Capillary length density, pc 
Capillary volume fraction, 
Total capillary number, Nc 
Capillary segment length, Lc 
Capillary blood velocity, Uc 
Capillary transit time, Tc 
Capillary oxygen pressure, po2 
Capillary length per neuron, Pc/Pn 
Cerebral blood flow, CBF 
Blood flow per neuron, CBF/p„ 
Oxygen consumption rate, CMR02 
Oxygen use per neuron, CMRo2/Pn 



171712 

Pc - V-^l' 

^ 1/2/3-7/24 
- -1/1/6+7/24 

P02 ~ 

Pc/Pn ~ 

CBF ~ 1/-^/6 
CBF/p„ ~ V 

CMR02 ~ v-^/^ 
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Table S2: Cerebral capillary and neural characteristics in mammals. 



Species 


2Rc 


fc 




Pc 




Pn 


Pel Pn 




\pm] 






[mm/ 


mm^] 


[lOVcm^] 


[//m] 


Mouse 


3.1 


0.015 


(b) 


880.0 




116.0 ^'^^ 


7.59 


Rat 


4.1 ('^^ 


0.014 




806.0 


ii) 


71.4 « 


11.29 


Cat 


5.1 ('^) 


0.021 


id) 


780.0 


ij) 


50.0 (™) 


15.60 


Dog 


4.5 










12.6 ('"^ 




Rhesus monkey 


5.6 (^) 


0.009 


if) 


329.9 


if) 


37.4 


8.82 


Human 


6.4 


0.023 


i9) 


219.0 


i9) 


38.5 


5.69 



The ratio pd Pn is the average capillary length per neuron. 

References: [14]; ^ [51]; [32]; (<^) [40]; [79]; [80]; [33]; [34]; « 

[78]; W [52]; « [24]; M [82]; ^ [41]. 
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Table S3: Arterial partial oxygen pressure and average cortical CBF. 



opecies 


arterial po2 


mm Hg 


'^J-'^ cortex 


PRP //I 
"^^^ cortex/ yn 


iviouse 


93.0 


iO / .4 ^ ' 


i.oD 


Rat 


93.0 


131.8 


1.84 


Rabbit 


75.0 


70.2 ('^) 




Cynomolgus 








monkey 


102.4 ('^) 


57.9 


1.51 


Dog 


107.0 






Rhesus monkey 


92.4 (-^^ 


52.5 


1.40 


Pig 


104.0 


59.5 




Human 


94.0 


42.6 (^'^ 


1.11 



The ratio CBF cortex /Pn is the average cortical blood flow per neuron (10 ^ mL/min). 

References: [70]; ^ [71]; [72]; ('^) [73]; [84]; [74]; ^s) [75]; [16]; W [83]; (^'^ 
[76]. 
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Appendix SI 

In this Appendix we want to show a formal step connecting Eq. (6) and Eq. (8) in 
the main text. We want to find the exponent 7, which satisfies the following equation: 
{\n{RjR;)f'^ = {RjRoV. 

Taking the natural logarithm of both sides leads to (2/3) ln(ln(i?c/-Ro)) = 7 ^^{Rc/ Ro)-i 
which yields Eq. (7) in the main text. 
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